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Abstract 

We show that within classical statistical mechanics without taking the thermodynamic limit, the 
most general Boltzmann factor for the canonical ensemble is a (/-exponential function. The only 
assumption here is that microcanonical distributions have to be separable from of the total system 
energy, which is the prerequisite for any sensible measurement. We derive that all separable 
distributions are parametrized by a mathematical separation constant Q which can be related to 
the non-extensivity g-parameter in Tsallis distributions. We further demonstrate that nature fixes 
the separation constant Q to 1 for large dimensionality of Gibbs T-phase space. Our results will 
be relevant for systems with a low-dimensional T-space, for example nanosystems, comprised of a 
small number of particles or for systems with a dimensionally collapsed phase space, which might 
be the case for a large class of complex systems. 
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INTRODUCTION 



There is a certain need in the physical, chemical, biological, social and economical sciences 
to understand the origin and ubiquity of power law-distributions. Many of these distributions 
appear to be g-exponentials upon closer inspection, which are defined as 

e% = [l + (l-q)x]M-4 , (1) 

for l+(l—q)x > 0. An appealing attempt for a general approach to the matter is to construct 
a thermostatistics, where with the use of classical principles, such as the variational principle 
and the second fundamental theorem of thermodynamics, one would be able to naturally 
derive these distributions. The canonical entropy associated to g-exponential distribution 
functions is the so-called Tsallis entropy 

where p is the the normalized energy density and dT indicates phase space integration. 
Tsallis entropy which is a generalization of Boltzmann-Gibbs (BG) entropy is in principle 
non-additive and non-extensive. Classical BG entropy is recovered in the limit q — > 1. 

In contrast to the vast amount of papers in non-extensive statistical physics, there has 
been relatively little effort to derive power-law distributions, in particular g-exponential 
distributions, from first statistical mechanics principles. To mention some fruitful work 
in this direction, in [3( it was shown that if phase space volumes are not covered in a 
Poissonian manner, the resulting entropy functional can be of Tsallis type. Within the 
framework of superstatistics 4] the inverse temperatures /3 are considered to fluctuate, such 
that a generalized Boltzmann factor, B = j df3f{(5)e~^ E arises. If / is a x 2 distribution, the 
natural distribution functions are g-exponentials. An other approach was taken in |5j], where 
the non-uniqueness of counting rules are discussed. It is shown that a propper modification 
of the log- counting rule yields Tsallis entropy instead of the usual BG entropy. In a recent 
paper fjj, we could show that a more general view is possible, in the sense that the most 
general Boltzmann factor that can be derived for the canonical ensemble is exactly a q- 
exponential. This result might be of relevance for systems with a low dimensionality of V 
phase space (Gibbs). Examples for such systems are nanosystems, which have gained recent 
interest due to being at the edge of technical accessibility, and for systems, whose phase 
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space has effectively collapsed in dimensions (like a fractal), due to long range interactions, 
aging, synchronization, as is the case in many complex systems. 



DERIVATION OF THE GENERAL BOLTZMANN FACTOR 

The following argument is solely based on the observation that any thermodynamic sys- 
tem which can be measured in equilibrium must be energy separable, i.e. thermodynamic 
quantities of the measured system must not explicitly depend on the energy of the total 
system, E. We consider a sample (observed system) in contact with a reservoir (bath). The 
energy of the sample is E\, the energy of the reservoir is E 2) such that the total system has 
a constant total energy E = E\ + E 2 . The number of microstates are u\(Ei) and co 2 (E 2 ) 
for the sample and the bath, respectively. The energy of the sample fluctuates around its 
equilibrium (extremal) value denoted by E*. Thermal contact of the two systems means 
that the Hamiltonian of the total system is H — Hi + H 2 and the partition function Z(E) 
is the convolution of the two microcanonical densities 

E 

Z{E) = J dEiUi(Ei)u 2 (E - Ei) with u^E,) = J dY t 5(Hi - . (3) 
o 

Following usual reasoning, we pass from the microcanonical to the canonical description. Its 
density p is given (up to a constant multiplicative factor) by 

p{Ei)=ui{E 1 )oo 2 {E-E l )Z-\E) . (4) 

Note, that this description is entirely dictated by the equations of motion. Assuming the 
existence of a unique extremal configuration at some E\ = E* defined by 5p = 0, leads to 
the definition of inverse temperature 

^11 _ g 1 (r\ 

Ul \E 1 =E t U2 \E 2 =E-E t - P k B T ' K ) 

The usual definition of entropy Si = k ln(cjj) implies that the extremal configuration is found 
where S — Si + S 2 is extremal with its associated temperature as defined above. Under 
which circumstances can one factorize the dependence of p on the total energy E, i.e. which 
microcanonical distributions allow for a separation of E into a multiplicative factor? The 
standard way to motivate the appearance of the Boltzmann term in the canonical ensemble 
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is a consequence of this ^-separation 

uj 2 (E-Ex) = enp(\n(uj 2 (E - Ex))) 

« exp(\n(u 2 (E)) - ^ ln(u; 2 ) Ex) (6) 
^ u 2 (E)exp{-/3E 1 ) . 

The approximation in Eq. © is exact for u 2 (E — Ex) being an exponential in E. However, 
this is not the most general way of separation. 

Now, to find the most general separation, we generalize the log function in Eq. © to 
some real function /, being twice different iable and with a well defined inverse The 
idea is to write u(E — Ex) = f^ 1 o f ou((E — E*) — (Ex — E*)) and to expand f ou around 
E — E*. If the energy E is separable from the system into a factor, then there must exist 
two functions g (factor) and h (general Boltzmann term) such that 
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(E-Ex) = g(u(E-E*)) h(x) , (7) 



with x := P(Ex — E*); to simplify notation we write u> := u(E — E*). Now use / to find the 
unknown functions g and h by expanding / o u to first order 

f(u(E-Ex))=f(g(0)h(x))~f(o)-Oxf'(u) , (8) 

which is justified for small x, i.e., the system being at or near equilibrium. The most general 
solution to this separation Ansatz is given by the family of equations (/, g, H)q, parametrized 
by a separation constant Q, and C and C 2 being real constants 

f(uj) = Cu^ + Ca 

g(oj) = oo (9) 
h(x) = [l-(l-Q)x] T ^ 

To prove this, set x = and h = h(0), so that Eq. (jHJ) yields f(g(uj)ho) = f(u>), which means 
g(uj) = ^. Without loss of generality set h = 1, and arrive at f(u)h(x)) = f(u>) — u)xf r (u). 
Form partial derivatives of this expression with respect to x and u>, and eliminate the f'(uh) 
term from the two resulting equations 

f'(u)h)h> = -f'(u) 

f'{ujh)h = (l-x)f'-uxf» 
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to arrive at the separation equation 

where Q is a mathematically necessary separation constant. The differential equation f — 
i (rr + l) = Q is straight forwardly solved to give the general Boltzmann term, h(x) = [l — 
(f — Q)x\ 1_Q , using h(0) = 1 to fix the integration constant. Equation —lujj^j = Q means, 
f{uj) = Cij3q^ 1-q + C 2 , with G\ and C 2 integration constants. / is strictly monotonous 
except for Q — 1, where it is constant. The term of interest in the canonical distribution 
can now be written as the generalized Boltzmann factor 

u 2 (E-E 1 ) =LU 2 (E-E*)[l-(l-Q)f3(E 1 -E*)] T ^ . (12) 

The usual Boltzmann factor Eq. © is recovered as the special case in the limit Q — > 1. If 
u! 2 is of power-form, tu 2 oc E 1 / 1 ^®, as is the case for a huge class of physical systems, Eq. 
(jT21) holds exactly. 

The separation constant Q is not specified at this level. What fixes Q7 The choice of 
thephysical system (Hamiltonian and characteristics of phase space) does. As an example, 
in pi we have shown that for a iV-particle Hamiltonian with pair-potentials, in D space 
dimensions 

N 2 N 



the following relation, fixing the separation constant, holds exactly, 

1 (a + 2)n 



i-Q 2a - 1 • < 14 » 

This equation establishes the connection between the interaction term in the Hamiltonian a, 
the dimensionality of phase space, n = DN, and the separation constant Q. From Eq. (|14|) 
it is immediately clear that for large systems the separation constant is always Q — > 1, i.e. 
the classical Boltzmann term © is recovered. For small systems, either due to low particle 
numbers or due to an effectively collapsed phase space dimensionality, a non-trivial Q / 1 
is expected. 



ON THE DEFINITION OF TEMPERATURE IN LOW-DIMENSIONAL SYSTEMS 

For large systems, the notion of equilibrium is well defined. However, for systems where 
we would expect Q ^ 1 the definition of temperature is not necessarily unique. In Eq. (J3J) 
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the definition of (3 or the inverse temperature was based on the extremal value of p, but this 
is not what is going on in a measurement. A measurement of temperature yields an expected 
value (due to averages of kinetic energy taken in the measurement process) and not the 
extremal value, i.e. the most likely one. To construct a theory which is consistent with the 
measured temperature of the sample system, and by - at the same time - keep the extremal 
principle, one can now ask to modify the definition of p — > p, such that the extremal value 
is obtained at the measured temperature. Effectively this amounts to an energy shift from 
equilibrium energy to expected energy, E* — > E. Such a modification could look like 

p(E 1 )=u 1 (E 1 )u) 2 (E-E 1 ) , (15) 

at Ei = Ei and where uJi(Ei) ~ E\ l . The idea is to identify E\ with the measured energy 
Ei = deep(e), and not with the equilibrium E* as before. The variation dp = leads to 
the relation 

¥ = > (16) 

A2 E — Ei 

which establishes the relation of the As, which are of course not independent. Now we note 
that 

Ei = [ 1 deeui(e)u 2 (E - e)Z(E)- 1 = Ef3{\i + 1, A 2 )//3(A 1; A 2 ) , (17) 
J 

with the usual beta functions, (3(a,b) = J Q dxx a {l — x) b . Substituting E\ into Eq. ([Tfi|) 
allows to compare the A exponents from the original system, Eq. with those of the 
energy shifted system, A, 

h = I as) 

A 2 /?(A 1 ,A 2 )//3(A 1 + 1,A 2 )-1 • 1 ; 

If we now assume that the exponent of the sample system is known (for example from a 
Hamiltonian), and unchanged in the shifted system, Ai = Ai, we observe that the shift from 
E* — > E, has as a consequence for the bath density, co , (£' 2 ) — > c<j(£ , 2 ) 2 , where z = ^- ^o^+i A2) • 
For the usual entropies this means 

Si = \n(ui) S 2 = ln{u z 2 ) (19) 

with the temperatures (ks = 1) 1/Tj = dSi/dE so that equi-temperature is assumed when 
both systems contain the expected energy, i.e. T\ = T 2 <->• SS = 0. Here we obtain 
non-additivity of entropies as a result of harmonizing the concept of temperature for small 
systems. 



In order to see that this gives consistent results we compute the expected temperature 
of both systems (Ti) and (T 2 ) in both, the most-probable state temperature interpretation 
and in the expected-state temperature interpretation. Let us think of the sample system 
(1) as the thermometer then the above definition implies that our thermometer shows the 
same average temperature obtained by multiple measurements whether we use the maximal 
or expected energy interpretation and (Ti) = (E/Xi)P(X 1 + 1, A 2 )//3(Ai, A 2 ) since we have 
chosen to modify only the exponent of the bath system. Now the bath shows the expectation 
of temperature, (T 2 ) = (E/\ 2 w)f3(\ 1: A 2 +l)//?(Ai, A 2 ) with w = 1 for the most-probable and 
w = z for the expected-state temperature interpretation. Consequently, the expected-state 
temperature interpretation is consistent with the BG expectation, 

<7\> Ai* /?(Ai + l,A 2 ) 

(T 2 > A 2 flAi.Aa + l) • 1 J 

Note that the expected-temperature definition involves knowledge about sample and bath 
system. On the other hand the most-probable-state temperature interpretation gives 

(Ti) _ A 2 /?(Ai + l,A 2 ) 

(T 2 ) A 1 /?(A 1 ,A 2 + 1) 1 J 

which in general will not equal unity. The corresponding notion of equilibrium allows to 
define the most-probable entropy merely upon the knowledge of the particular micro-systems 
separately. Since we only actually measure the sample system this temperature difference can 
not be observed. Averaging over multiple temperature measurements does not converge to an 
equailibrium temperature. The way of how to average multiple temperature measurements 
in order to obtain the equilibrium temperature in terms of the most-probable state therefore 
is dual to the manipulation of the bath entropy, constructed to match the expected-state 
interpretation. 

It becomes evident that the usage of expected and most probable energies can both be 
consistently used in order to define the macroscopic notion of equilibrium. Which of the 
two possibilities is realized experimentally is defined by the choice of the experimentator. 
Intuitively we would think that in terms of keeping records of repeated experiments the 
expected-state temperature interpretation seems more natural since we can directly identify 
averages on multiple measurements of the temperature with the equilibrium temperature 
itself. Whatever procedure we fix, we can not escape modifying either, the average on the 
replicas, i.e. multiple measurements, or the definition of the bath entropy. The asymmetric 
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definition of the entropy does not appear explicitly in the canonical variational principle, 
based on single particle entropies where the same mechanisms are controlled by constraints. 

CONCLUSION 

Based on very general assumptions we have derived that the most general Boltzmann 
factor in the canonical ensemble is a g-exponential. We have shown that this result might 
be relevant for systems of relatively small phase space, which might be realized for several 
types of complex systems. We comment on the temperature definition in small systems 
and discuss its consequences. Let us stress again, that all presented arguments are strictly 
based on Hamiltonians and on the variational principle, we never leave the field of classical 
statistical mechanics, except for not taking the thermodynamic limit. 
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